Positive Wigner functions render classical simulation of quantum computation efficient 
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We show that quantum circuits where the initial state and all the following quantum operations can be rep- 
resented by positive Wigner functions can be classically efficiently simulated. This is true both for continuous- 
variable as well as discrete variable systems in odd prime dimensions, two cases which will be treated on entirely 
the same footing. Noting the fact that Clifford and Gaussian operations preserve the positivity of the Wigner 
function, our result generalizes the Gottesman-Rnill theorem. Our algorithm provides a way of sampling from 
the output distribution of a computation or a simulation, including the efficient sampling from an approximate 
output distribution in case of sampling imperfections for initial states, gates, or measurements. In this sense, this 
work highlights the role of the positive Wigner function as separating classically efficiently simulatable systems 
from those that are potentially universal for quantum computing and simulation, and it emphasizes the role of 
negativity of the Wigner function as a computational resource. 



What renders a quantum computer a superior computa- 
tional device? Where is the precise boundary between clas- 
sically efficiently simulable problems and ones for which this 
is no longer possible? Despite a significant research effort 
and partial progress ifl UIoll . these questions are still largely 
open. Quantum correlations surely play a role in one way or 
the other in quantum computers and simulators outperforming 
their classical counterpart. For example, if the entanglement 
is - in a precise sense - too low in a pure-state computation 
with respect to any bi-partite split, then one can classically ef- 
ficiently simulate the dynamics ||2|-|6t]. In measurement-based 
computing specifically ll ill 1211 . where the resource character 
of entanglement is particularly manifest, states can be too little 
entangled 0], but also in a sense too entangled J3l. Possibly 
the singularly most important result on classical simulatability 
of quantum computers is the Gottesman-Knill theorem, stat- 
ing that stabilizer circuits consistingof Clifford gates only can 
be classically efficiently simulated I1 1 L I SL 1 1 3fl . Though, adding 
almost any further gate will render the Clifford gate set uni- 
versal for quantum computing. Similarly, Gaussian operations 
for continuous-variable systems can be efficiently simulated 
CS El- 
Here we present a generalization of the Gottesman-Knill 
theorem, stating that one can efficiently classically sample 
from quantum circuits starting from product states with a posi- 
tive Wigner function, applying quantum gates that have a pos- 
itive Wigner function (in a sense made precise below) and per- 
forming measurements associated with positive Wigner func- 
tions. This result holds true both for discrete variable systems 
where the constituents have odd prime dimension (it can eas- 
ily be generalized to arbitrary odd dimension) as well as for 
continuous-variable systems so common in quantum optics. 
In fact, these two situations can be treated on exactly the same 
footing - since at the root of the remarkably simple argument, 
we exploit the structure of the phase spaces, discrete as well 
as continuous. The relative elegance of the approach draws 
from the observation that the expressions for the discrete and 
continuous description are identical. 

In this sense, negativity of the Wigner function is identi- 



fied as a computational resource. Indeed, it is usually seen 
as an indicator of "non-classicality" Il8l420ll (implying the P- 
function to be non-positive). Here we see that it is exactly this 
negativity that is needed if a quantum computer or simulator is 
to outperform its classical counterpart. This results adds meat 
to the notion of such states being "classical", quite similar to 
the situation of states with positive Wigner functions and ho- 
modyne measurements being unable to violate a Bell inequal- 
ity l2lll or be useful in magic state distillation 112211 . We also 
comment on the converse direction, that quantum computing 
is possible with resources having Wigner functions arbitrarily 
close to being positive. While negativity is a necessary re- 
source, one does not need to have "a lot of it" and it can be 
used in a dilute form. 



DISCRETE AND CONTINUOUS PHASE SPACES 

We start by discussing the less commonly addressed phase 
space structure of finite-dimensional quantum systems (com- 



pare Refs. 1231 1241 126H ). and shift the emphasis to infinite- 
dimensional ones later. We assume the local dimension d to be 
an odd prime, merely to avoid technicalities required in non- 
prime dimension (but they can be treated on the same footing 
l26ll ). In this case the phase space of a single d-level systems 
is TLd x %a, so that it can be associated with a d x d cubic 
lattice. Given an orthonormal basis {|0), |1), . . . \d — 1)}, one 
can define the shift and boost operators as the generalizations 
of the familiar Pauli matrices by 

x(q)\x) = \x + q mod d), z(p)\x) = uj px \x), (1) 

where u> = e 27 "/ d is a <i-th root of unity, arithmetic being 
modulo d, and x = 0, . . . , d — 1. The fundamental tools of 
every quantum phase space representation are the so called 
Weyl operators. For discrete systems they are given by 



w(q,p) = uj 2 pq z(p)x{q)i 



(2) 



where 2 1 = (d+l)/2 is the multiplicative inverse of 2 (again 
modulo d). These operators form a group, the Heisenberg- 
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Weyl group, and are the main ingredient for representing 
quantum systems in phase space. 

Let us now turn to continuous systems, and focus on a 
single mode first, associated with canonical coordinates or 
"quadratures" of position Q and momentum P. The associ- 
ated phase space is now R 2 . Indeed, the continuous version 
of the Weyl operator can be analogously given in terms of the 
previous formulas (fl3 and ©, with the only difference us- 
ing standard arithmetic. In fact, with a different choice of the 
phase factor, namely u = e\ the Weyl operator in Eq. (01 
takes the form 



three important properties of the Wigner function: 



w(q,p) — exp(igP — ipX), 



(3) 



which is the familiar displacement operator well known in 
quantum optics. From now on we are going to use a unique 
notation defined in such a way to be consistent with discrete 
and continuous phase space representations. For this reason, 
in Table |T] we introduce a set of symbols which are valid in 
both settings. 



Symbol 


Discrete 


Continuous 




e 2 ™ /d mod 


d 


e* 


w(q,p) 


uj- 2 ~ lpq z(p)x 


(?) 


exp(iqP — ipX) 


r 


(qi,pi, . . -,q n ,Pn) 


ez 2 d n 


(qi,pi,...,q n ,p n ) € R 2 ™ 


£ 






J rSRaiZ o(r) ar 


c 


2d 




2ir 



TABLE I. Table of symbols providing a unified notation valid for 
both discrete and continuous-variable systems. The function size(r) 
gives the number of constitutents of r, so that the sum (or the integral) 
is correctly defined also for local subsystems. 



tv(p) 
tr(AB) 



W p (r) = 1, 

W A (r)W B (r), 



(5) 



(6) 



W AT (r) = W A (Ar), A = diag(l, -1, 1, -1). (7) 

By virtue of Hudson's theorem, the only pure states having 
a positive Wigner function are Gaussian states 12711 for con- 
tinuous systems 12811 and stabilizer states for odd local dimen- 
sion d [12611. Unitary operations preserving the Gaussian form 
of the Wigner function are Clifford operations in discrete sys- 
tems and Gaussian operations in continuous systems. Those 
operations admit, via the Jamiolkowski isomorphism, a de- 
scription in terms of positive Gaussian Wigner functions. 

For mixed states the situation is more complex. Surely, con- 
vex combinations of Gaussian and stabilizer states will have 
positive Wigner functions. But there also exist states with a 
positive Wigner function which cannot be represented in this 
form 1 26[ 231 ■ Analogously there are quantum operations ad- 
mitting a positive Wigner representation but which are not a 
convex combination of Clifford/Gaussian unitary operations. 
In this work we will focus on the simulation of quantum algo- 
rithms involving states, operations and measurements all de- 
scribed by positive Wigner functions. This corresponds to a 
more general scenario of which the Gaussian setting is a par- 
ticular case. In this sense our result can be viewed as an ex- 
tension of the Gottesman-Knill theorem. Quite remarkably, 
our method it is completely independent from any Gaussian- 
ity property. 



Operations having positive Wigner functions 



The Wigner function of an operator O acting on n discrete 
or continuous-variable systems is defined as 

W {r) = (c/2)""tr {w{r)Ii® n w{r)^ O) . (4) 

where w(r) = w{qi,p\) (g> • • • <g> w(q n ,p n ) and II is the single 
system parity operator: For discrete systems, on state vec- 
tors, this parity operator acts as II : \x) M- | — x mod d) 
for continuous-variable systems as (II^)(ar) = ip(—x). This 
function has the structure of W : H?^ — > R for discrete sys- 
tems (compare also Ref. dl) and W : R 2n R for contin- 
uous systems. 



Properties of the Wigner function 

The Wigner function of a state p is normalized and can be 
interpreted as the quantum analogue of a phase space distribu- 
tion with the peculiar property of being not necessarily posi- 
tive in all its domain. In the next sections we are going to use 



We have already introduced the Wigner function of a state 
or more generally of an operator. We still have to properly 
give a phase space description to operations and measure- 
ments. We now define the Choi matrix f of a completely 
positive map F. We define 

/ = (l®f)|w)(w| (8) 

where \uj) = (Y^j—i \ jij))® n f° r discrete systems, which is 
up to normalization the state vector of the maximally entan- 
gled state. For continuous systems, if the limit exists and is of 
trace class, we set 

/= ®F)\lo s ){lo s \ (9) 

where \uj s ) = (Y^jLi s3 j))^ n - The Choi matrix is posi- 
tive, as a consequence of F being completely positive, and it 
is supported on a Hilbert space having the natural structure of 
H = H ln T-l° ut . We denote the partially transposed matrix 
with respect to subsystem "out" with the symbol / r . We will 
say that the completely positive map F "has a positive Wigner 
function" if the Choi matrix / has a positive Wigner function 
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Wf. It is easy to see that if / has a positive Wigner function, 
then the same is true for / r as well. The application of an op- 
eration to a state p out = F(p m ) is related to the Choi matrix 
via a partial transposition T and a partial trace 

p° ut = tr in ((p in ® l)/ r ). (10) 

In phase space this is reflected by 

W(r out ) = c 2n j> W r (r out ;r in )W(r in ), (11) 

where r ln , r out £ E 2 ™. Trace-preserving operations are such 
tr ou t (/) = 1, i-e. 

c 2n <f W f r(r out ;r ia ) = l, (12) 

J r out 

for all r ln . This means that, if the function is positive, c 2n Wjr 
can be interpreted as a classical stochastic matrix. This will be 
a key property for the classical sampling algorithm. For dis- 
crete systems, the Wigner function associated with the identity 
operation is given by the Kronecker delta 

W ]u)){u] r =c- 2n 6(r°» t ;r in ). (13) 

For continuous systems, one can also consider Wigner func- 
tions associated with operations for which the Choi matrix is 
not of trace class, such as when F is the identity operation. In 
this case, we have in the sense of distributions 

W lUs)(Us] r^c- 2n 5(r out ;r in ). (14) 

An important subcase of operations having a positive 
Wigner functions is constituted by Clifford/Gaussian unitaries 
Us.d- They map Weyl operators onto Weyl operators under 
conjugation, so that Us,d,w(r)Ug d = w(Sr + d) is again a 
valid Weyl operator. Here, S 6 Sp{2n, 7L d ) and v e 7L 2 d n 
for discrete systems while S <E Sp(2n, R) and v e R 2n for 
continuous systems. From Eq. ( TBl i. one can deduce that the 
associated Choi matrix has the following Wigner function rep- 
resentation for discrete systems l26ll 

W f v sd = c- 2n S{Sr out + d;r in ), (15) 

and similarly for continuous systems in the sense of distri- 
butions. Physically, Clifford and Gaussian operations are 
of utmost importance in the laboratory since they can often 
be realized with simple experimental settings. Since Clif- 
ford/Gaussian operations and convex combinations have a 
positive Wigner function, our result covers also this setting. 

Measurements with positive Wigner functions 

We finally have to give a phase space description of a gen- 
eral measurement acting on single tensor factors. Given a gen- 
eral measurement defined by a set of positive operator valued 



measures (POVM) {M k } satisfying £)fc=i M k = 1, we as- 
sociate to each operator M k the respective Wigner function 
W Mk : 7!? d ->■ R for discrete systems and W Mk : R 2 R 
for continuous systems. Given a state p the probability of 
getting the measurement outcome k on subsystem I will be 
P{k) = tr{pMf, ) and hence, for a single system, 

P{k) = c 2 j>W p {r)W Mk {ri). (16) 
Moreover, we also have 

K 

c 2 J2 Wm ^ = 1 (17) 

fe=i 

for all r both in the discrete and the continuous setting, which 
means that, if the functions c 2 WM k are positive, then they can 
be interpreted as probabilities in the variable k for each value 
of r. This fact will be also crucial in the simulation algorithm. 

STATEMENT OF THE PROBLEM 

Having laid out the formalism to be employed, we are now 
in the position to precisely state the problem at hand: We al- 
low for general quantum circuits of the following form, again 
both in the continuous or discrete setting. 

• We consider initial product states p = pi <g> • • • <g> p n . 

• To this initial state, a sequence of trace preserving quan- 
tum channels Ft o . . . F-x o F± is applied (not necessarily 
unitary), each of them supported on at most m subsys- 
tems at a time. 

• Finally, local measurements are performed on each in- 
dividual subsystem defined by some local projective 
positive operator valued measures {M k }. 

Here n denotes the number of constituents, while t is the 
depth of the circuit, so the number of gates or local operations 
applied. The classical simulation should scale polynomially 
with respect to these two parameters. We remark that initial 
states can be mixed and the applied gates will in general be 
non-unitary. This is probably the most general model of a 
non-adaptive quantum algorithm, i.e. where the sequence of 
gates and measurements is fixed. Non-adaptivity has be cho- 
sen just for simplicity of the exposition but the algorithm can 
be easily extended to the adaptive case. 

A run of the quantum algorithm will provide one list of out- 
comes fci, &2, . . . fc„, one for each measured subsystem. The 
probability of a given list of outcomes to occur is 

P(h, . . . , kn) = tr ((M kl ® M k2 . . . M kn )(F t o...F 2 o F x ){p)) . 

(18) 

Using quantum systems, one can hence sample from the distri- 
bution P. Classically, the problem to be solved is again a sam- 
pling problem: The quantum circuit can be classically simu- 
lated if there is a classical algorithm that is efficient in t and n 
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and provides, in each run, a list of outcomes k\ , . . . , k n drawn 
from (approximately) the same probabilities of the quantum 
circuit given in Eq. ( fT8l ). Note that we do not require to really 
compute the probabilities of all the possible outcomes (simu- 
lation in a stronger sense), but we just want a classical algo- 
rithm to be efficient in sampling from the distribution defined 
by the quantum circuit, in the sense that output strings of the 
classical and quantum machines are drawn from the same (or 
approximately the same) probability distribution. This weaker 
sense of simulation is anyway enough to exclude any possible 
speedup of the quantum algorithm with respect to the classical 
one, and hence identifies negativity of the Wigner function as 
a necessary resource in quantum computing. We will also take 
into account possible errors in the classical sampling, show- 
ing that those errors cannot strongly affect the final probability 
distribution. 

Phase space representation 

First note that the Wigner function of the input product state 
p = pi ® ■ • • ® p n is given by a product function W p with 

W p (r) = W p ,( ri )W P2 (r 2 )...W Pn (r n ), (19) 

where W Pl is associated with the subsystem pi. To each gate 
F t can be associated the Wigner function Wtr . The Wigner 
functions of the POVM associated with the outcome k and 
performed on subsystem I are denoted as Wn k with local 
phase space coordinates r; . By sequentially applying Eq. ©, 
we can express the outcome probabilities of the quantum cir- 
cuit as 

P(fc 1 ,...,fc„) = c 2n(t+1) / (20) 

J r (t) )r (t-l) ... |r (0) 

xW Mki (r^)...W Mkn (rW) 

x W f r (r (t) ; r {t ~ 1] )... W f r (r (1) ; r (0) )W P (A Q) ) , 

where the subscripts in the coordinates indicate subsystems 
while superscripts indicate the integer time steps associated 
with the sequential application of the gates. This is just a for- 
mal phase space description the quantum circuit completely 
equivalent to the operator representation given in Eq. (TT~8l >. 
Up to now no assumption has been made on the the type of 
resources and operations. 

EFFICIENT CLASSICAL SIMULATION 

We now make the following two assumptions: 

• The Wigner functions of the input state, the gates, and 
the POVMs are positive, 

• it is possible to draw phase space points according to lo- 
cal probability distributions associated with local states, 
local gates and local measurements. 



Later, the second of this assumptions will be relaxed allowing 
for eventual classical sampling errors. 

Observation 1 (Efficient classical simulation of circuits) 

For any n and t one can sample classically from the 
distribution P in poly(n, t) time. 

This is done as follows. 

• Step 0: Draw a phase space point f^> according to the 
input Wigner function W p (.), 

• Steps j = 1, ...,£: draw a phase space point 
f( J ) according to the probability distribution 

c 2n W f v(.-f^-^). 

J 3 

• Step t + 1: Finally, draw a measurement outcome 
k\ , k%, . . . , k n according to the probability distribution 

P(h, ...,k n )= c 2n W Mkl (f[ t] ) • . . W Mhn (fW). (21) 

Proof: This algorithm is simply a classical stochastic pro- 
cess and directly from the law of conditional probabilities we 
have that, in the final step t+1, the probability of getting the 
outcome ki, k 2 , ■ ■ • , k„ is given by Eq. d20l i. So one run of 
the classical algorithm is completely equivalent to one run of 
the quantum circuit. Moreover, the efficiency of the classical 
procedure with respect to n and t follows from the following 
observations: 

• Step is efficient because, since the initial state is a 
product, we merely have to draw n independent sub- 
system phase space vectors. 

• Steps 1,2, ... ,t are efficient because, for each gate in 
step j, we draw a phase space vector associated with at 
most m local subsystems and leave the complementary 
coordinates invariant. This follows from the structure 
of Wtr factorizing as (as a Kronecker delta or in the 
sense of distributions) 

c^rfrV'- 1 ) = c 2m w%°^\riy- l )8^ c yj 1 ), 

3 J j 

(22) 

where tl are the local coordinates of the m subsystems 
involved in the gate, while rc are the complementary 
ones. For every given input vector r 3_1 one has to 
draw a vector r J L with respect to the local distribution 
c 2m Wf° and just leave the other coordinates invari- 

3 

ant r^T 1 h-> r J G . 

• Step t + 1 is efficient because, since the final measure- 
ments are local, we draw n independent outcomes asso- 
ciated to each subsystem. 

As a final remark we observe that the positivity of all the 
Wigner functions and the properties given in Eqs. (|5), (fT2l 
and ( TPTI i are crucial for all the functions appearing in the clas- 
sical algorithm to be interpreted as probability distributions. 
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If this is not the case, the classical algorithm cannot be ap- 
plied. However, using conditional probabilities, the algorithm 
can easily accommodate adaptive later steps based on earlier 
measurement outcomes. 



ROBUSTNESS TO SAMPLING ERRORS 

In this section we discuss the robustness of the previous 
method with respect to possible errors in the classical sam- 
pling from positive Wigner functions. Suppose that in the 
classical algorithm one is able to efficiently sample phase 
space vectors only from imperfect probability distributions 
which are close to the ideal ones. These imperfect Wigner 
functions will correspond to imperfect density operators p' , 
gates F' with Jamiolkowski isomorphs /' and POVM ele- 
ments M' k . We will show that the classical algorithm is ro- 
bust with respect to accumulated sampling errors. What is 
more, these errors can be naturally linked to errors in the 
physical quantum system. For brevity, we focus on the dis- 
crete case here, while the continuous situation can be treated 
on the same footing. For simplicity of notation, in turn, we 
denote the stochastic matrix associated with the k-th step by 
Q{k) _ c 2n W f r(.; .) and its (stochastic) approximation by 

•> k 

Q 1 ^ ■ || - ||ooj ||-||i denote the usual respective matrix p-norms. 

Observation 2 (Efficient approximation of circuits) For 

any n, t, and any e > 0, consider a sampling from the 
distribution P' obtained from initial states, operations, and 
measurements deviating from those in Observation 1 in that 

\W Pl (n)-W p[ ( n )\ < e, \W Mkl {n)-W M - ki (n)\ < e, (23) 

for I = 1, . . . , n and all r;, and \\Q^ k ' — Q'^\\oo < sfor k = 
1, . . . ,t. Then one can sample from P' with \\P — P'||oo < 
s poly(n, t) in poly(n, t) time. 

Note that the above estimates are also valid if one has trace- 
norm bounds for all states as well as operator norm bounds for 
the POVM elements. Using the duality of the operator norm 
and the trace norm, one finds for all r/ 



IWM-wMiKd-^n-p'th 



(24) 



for I = 1, ... ,n. Similarly, from the definition of the parity 
operator, we get 



\W Mki {n)-W M . k {n)\ < \\M kl -M' kl \\ c 



(25) 



for every I = 1, . . . , n, any r;, and each outcome. 

Proof: Iteratively inserting and subtracting terms and using 
the triangle inequality as well as the fact that \W PI (r;)| < 1 
for all states pi and all r/, we get for all r 

n 

\W p (r) - W p ,(r)\ < £ \W pl (n) - W p[ {n)\ < ne. (26) 
i=i 



Using the triangle inequality several times, one finds 
i \ /i 



sup 



k=t 



\ Q (k) w p M - II Q w p' M 



< 



(*0 



fe=t 



k=t 



k=t 



sup 



(27) 



\Q' (k) {W p -W p .) (r) 



k=t 



Using the sub-multiplicativity of the operator norm, 



\y[q 

k=t 



< 



n 

fc=i 



(28) 



as well as the fact that for all states p and all r we have that 
\W p (r)\ < 1, we find 



sup [Y[Q {k} W p (r)- (Y[Q' {k) W p ,)(r) 

r ^fc=i ' ^k=t 

t 

<Y \\Q {k) -Q' [k) +ne<(t + n)e. 

II oo 

Using that 



(29) 



1=1 



1=1 



(30) 

and combining these results, one arrives at the above observa- 
tion. 



SUMMARY AND OUTLOOK 

At this point, one may ask "how negative the Wigner func- 
tion has to be" in order to allow for universal quantum com- 
puting. Surely, one can perform universal quantum computing 
with positive Wigner functions initial product states, as well as 
with circuits of quantum gates and local measurements, all of 
which equipped with a Wigner function being arbitrarily close 
to being positive; so for each element of a circuit, the answer 
is "not much": For every e > and d = 3, there exist fami- 
lies of initial states p\ , . . . , p n with a positive Wigner function 
each, families of local gates with isomorphs /i, ■ • • , /t sup- 
ported on a number of sites constant in n and local POVMs 
Mfc x , . . . , Mfc re such that 



W fj (r) > -e, W Ml {n) > 



(31) 



for I = 1, . . . , n, j — 1, . . . , t and all ri,r, giving rise to 
circuits universal for quantum computing. This observation 
is obvious from the fact that product stabilizer measurements 
have a positive Wigner function. What is more, unitary gates 
can be "diluted" such that each gate is close to the identity 
operation. Also, by encoding the output, POVMs arbitrarily 
close to the identity can be used. Similar schemes can also de- 
veloped based on encoded cluster states as in Refs. Il2ll3lll . 
Note that this is not in contradiction with Observation 2, since 
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despite the error growing slowly, it grows beyond all bounds 
in an unbounded computation using imperfect gates. A les- 
son to be learned from these rather obvious examples is that it 
appears to be a very fruitful enterprise to meaningfully quan- 
tify the negativity as a resource in terms of a proper resource 
theory, as it has happened in entanglement theory ll32ll or in 
statistical physics A3 311 . 

In this work, we have shown that the negativity of the 
Wigner function can be grasped as a resource in quantum 
computing and simulation: If the basic elements of a circuit 
exhibit a positive Wigner function, the probability distribution 
of the quantum computation can be efficiently sampled. This 
remains true if one can only approximately implement each 
gate, in that the errors made scale favorably. In this sense, our 
result generalizes the Gottesman-Knill theorem for sampling 
outcomes of circuits, in a way where continuous and discrete 
systems are treated on exactly the same footing. We hope that 
the present approach stimulates further work on identifying 
the boundary between classically efficiently simulable quan- 
tum systems and those universal for quantum computing. 
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APPENDIX 

Preparation of random variables with given Wigner functions 

Here we briefly sketch how to draw phase space points ac- 
cording to given Wigner functions, both for discrete as well as 
for continuous-variable systems. There are several standard 
methods of how to proceed here, e.g., the rejection method or 
the inverse transform sampling method. We will be specific 
for sampling in R s , s G N, with a similar reasoning hold- 
ing for W d as well. Denote for a given non-negative Wigner 



function W D : R s 



R with 



F : R s 



R 



(32) 
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the cumulative distribution function, r G R 2s . This can be Then, given 2c? uniformly distributed random variables with 
written as realizations ui, . . . , U2 S > take 

n = (34) 

r 2 = F R2 ) Ri (u 2 ), (35) 
F(n, . . . , r 2s ) = F Rl (n)F R2lRl (rain) (33) r2s = F *l ^ ' ' • ^-iK,). (36) 



x 



-P 1 fl 2a |-Ri...-R2 i ,-i ( r 2s| r i ■ ■ • r 2s-i)- This will then be distributed according to W p 



